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ABSTRACT 

O"' 

o : 

, If Gamma-Ray-Bursts (GRBs) occur at high redshifts, then their bright afterglow 

emission can be used to probe the ionization and metal enrichment histories of the 
'-^ | intervening intergalactic medium during the epoch of reionization. In contrast to 

pL< ' other sources, such as galaxies or quasars, which fade rapidly with increasing redshift, 

{Nj ■ the observed infrared flux from a GRB afterglow at a fixed observed age is only a 

weak function of its redshift. This results from a combination of the spectral slope 
of GRB afterglows and the time-stretching of their evolution in the observer's frame. 
Assuming that the GRB rate is proportional to the star formation rate and that the 
characteristic energy output of GRBs is ~ 10 52 ergs, we predict that there are always 
~ 15 GRBs from redshifts z ^ 5 across the sky which are brighter than ~ 100 nJy at 
<^ ' an observed wavelength of ~ 2^m. The infrared spectrum of these sources could be 

taken with the future Next Generation Space Telescope, as a follow-up on their early 

■ X-ray localization with the Swift satellite. 

Q-f 

6 ■ 

Subject headings: gamma rays: bursts — ISM 

>'■ 

■ 1. Introduction 

The past decade has been marked by major observational breakthroughs concerning the 
properties of the Gamma Ray Burst (GRB) sources. The Burst and Transient Source Experiment 
(BATSE) on board the Compton Gamma Ray Observatory (Meegan et al. 1992) showed that the 
GRB population is distributed isotropically across the sky, and that there is a deficiency of faint 
GRBs relative to a Euclidean distribution. These were the first observational clues indicating 
a cosmological distance scale for GRBs. The localization of GRBs by X-ray observations with 
the BeppoSAX satellite (Costa et al. 1997) allowed the detection of afterglow emission at optical 
(e.g., van Paradijs et al. 1997) and radio (e.g., Frail et al. 1997; Frail, Waxman &: Kulkarni 1999) 
wavelengths up to more than a year following the events (Fruchter et al. 1999; Frail, Waxman & 
Kulkarni 1999). The afterglow emission is characterized by a broken power-law spectrum with a 
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peak frequency that declines with time. The radiation can be well modeled as synchrotron emission 
from a decelerating blast wave, created by the GRB explosion in an ambient medium, plausibly 
the interstellar medium of galaxies (Waxman 1997; Wijers & Galama 1999; Meszaros 1999; but 
see also Chevalier &; Li 1999). The detection of spectral features, such as metal absorption lines 
(Metzger et al. 1997), in some optical afterglows allowed an unambiguous identification of the 
cosmological distance scale. By now, the redshift of almost a dozen GRBs has been identified 
either through the detection of absorption features in the afterglow spectra or of emission lines 
from host galaxies (see Fig. 3 in Kulkarni et al. 2000). 

The central engine of GRBs is still unknown. Since the inferred energy release is comparable 
to or higher than that in supernovae, most popular models relate GRBs to stellar remnants, 
such as neutron stars or black holes (see, e.g., Eichler et al. 1989; Paczyhski 1991; Usov 1992; 
Mochkovitch et al. 1993; Paczyhski 1998; MacFadyen & Woosley 1999). Recently it has been 
claimed that the late evolution of some rapidly declining optical afterglows shows a component 
which is possibly associated with supernova emission (e.g., Bloom et al. 1999; Reichart 1999). If 
the supernova association will be confirmed by detailed spectra of future afterglows, the GRB 
phenomenon will be linked to the terminal evolution of massive stars. 

Any association of GRBs with the formation of single compact stars implies that the GRB 
rate should trace the star formation history of the universe (Totani 1997; Sahu et al. 1997; Wijers 
et al. 1998; but see Krumholz, Thorsett &: Harrison 1998). Owing to their high brightness, 
GRB afterglows might be detected out to exceedingly high redshifts. Similarly to quasars, 
the broad band emission of GRB afterglows can be used to probe the absorption properties of 
the intergalactic medium (IGM) out to the epoch when it was reionized at a redshift z ~ 10 
(Miralda-Escude & Rees 1998; Loeb 1999). Lamb & Reichart (1999) have extrapolated the 
observed gamma-ray and afterglow spectra of known GRBs to high redshifts and emphasized 
the important role that their detection might play in probing the IGM. Simple scaling of the 
long-wavelength spectra and temporal evolution of afterglows with redshift implies that at a 
fixed time lag after the GRB in the observer's frame, there is only mild change in the observed 
flux at infrared or radio wavelengths with increasing redshift. This results in part from the fact 
that afterglows are brighter at earlier times, and that a given observed time refers to an earlier 
intrinsic time in the source frame as the source redshift increases. The mild dependence of the 
long-wavelength flux on redshift is in contrast with other high-redshift sources such as galaxies or 
quasars, which fade rapidly with increasing redshift (Haiman & Loeb 1998; 1999). The "apparent 
brightening" of GRB afterglows with redshift could be further enhanced by the expected increase 
in the mean density of the interstellar medium of galaxies at increasing redshifts (Wood & Loeb 
1999). 

It therefore appears natural to use GRBs as an important tool in probing the high-redshift 
universe and its star formation history (Blain & Natarajan 1999). Since GRBs are rare, all-sky 
searches for their early 7-ray emission are needed before follow-up observations at much longer 
wavelength are conducted. 
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In this paper we model the emission from GRB afterglows and follow their number counts at 
different wavelengths as a function of redshift. For simplicity, we assume that the ambient gas 
surrounding GRB sources is the interstellar medium of their host galaxies. In §2 we review our 
model for both the relativistic and non-relativistic stages of the afterglow emission, as well as a 
simple prescription for the redshift evolution of the interstellar medium of galaxies. The numerical 
results and their dependence on model assumptions are discussed in §3. Finally, §4 summarizes 
our main conclusions and describes the implications of our results. 



2. Afterglow Emission 

2.1. Relativistic Regime 

GRB afterglows can be modeled as synchrotron emission from a decelerating relativistic blast 
wave, created by the GRB explosion in an external medium. For a point explosion in a uniform 
medium, the shock structure in the highly-relativistic regime is described by the Blandford 
& McKee (1976) self-similar solution. The synchrotron emission originates from a power-law 
population of shock-accelerated electrons in a strong magnetic field (Waxman 1997; Meszaros, 
Rees, & Wijers 1998; Sari, Piran & Narayan 1998), with both the electrons and the magnetic field 
having close to equipartition energy densities (see Medvedev & Loeb 1999 for the possible origin of 
the magnetic field). For simplicity, we focus our discussion on a spherically-symmetric explosion in 
a uniform ambient medium, assumed to be the interstellar medium of the host galaxy, although at 
least some GRB blast waves might be expanding into the stellar wind of the progenitor star (Dai 
& Lu 1998; Meszaros, Rees &; Wijers 1998; Chevalier &: Li 1999). The dependence of our results 
on the ambient medium density will be discussed in § 3.3. We assume that the shock-accelerated 
electrons have a power-law distribution of Lorentz factor, 7 e , with a minimum Lorentz factor j m 
(see, e.g., Sari, Piran & Narayan 1998). We also define 7 C as the threshold Lorentz factor below 
which electrons do not lose a significant fraction of their energy to radiation. In the fast cooling 
regime, when 7 m > j c , all of the electrons cool rapidly down to a Lorentz factor ~ 7 C and the flux 
observed at the frequency v is given by, 

Fv = F Vm l {u/uc)- 1 ' 2 v m >v>v c , (1) 

where F Um is the observed peak flux, v c = v{j c ) and v m = v{j m ) are the characteristic synchrotron 
frequency calculated at 7 C and j m respectively and p is the power-law index of the electron energy 
distribution. A typical value of p ~ 2.5 often fits both the GRB and the afterglow observations 
(Sari, Piran & Narayan 1998; Kumar & Piran 1999). 

When 7 C > 7 m , only electrons with 7 e > 7 C cool efficiently. In this slow cooling regime the 
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observed flux is, 

!W^m) 1/3 v m > V, 

(v/v m )- {p - 1)/2 V C >V> U m , (2) 

(Vc/Vm)- (P - 1)/2 (V/Vc)- P/2 V>V C . 

Typically, synchrotron self-absorption results in an additional break, at a frequency ^ 5GHz. We 
will not consider the low-frequency regime below this break in our discussion. Assuming a fully 
adiabatic shock (Sari, Piran & Narayan 1998), 

v c = 2.7 x W 12 e B 3/2 E- 2 1/2 niH d 1/2 (l + z)^ 2 Hz, (3) 

v m = 5.7 x 10M 2 ^V /2 (1 + z ) 1/2 Hz > ( 4 ) 

F„ m = 1.1 x 10 5 e^ /2 E 52 nJ /2 cZ 2 - 8 2 (l + *) ^Jy. (5) 

The redshift dependence results from the fact that the radiation emitted by a source at a redshift 
z at the frequency v s over a time Ai s , will be observed at z = at a frequency f D = za,/ (1 + z) over 
a time At Q = (1 + z)At s . Here eg and e e are the fraction of the shock energy that is converted to 
magnetic fields and accelerated electrons, respectively, for which we adopt the values €b = 0.1 and 
e e = 0.2 (see, e.g., Waxman 1997); E = .E52IO 52 erg is the energy of the spherical shock; n = n\ 
1 cm -3 is the mean number density of the ambient gas; t = t& \ day is the time lag since the 
GRB trigger, as measured in the observer frame, and di = c?28l0 28 cm is the cosmology-dependent 
luminosity distance. For a flat universe (Qq = VL m + 17a = 1), the luminosity distance can be 
written as: 

dz', (6) 



d L = (l + z) j\l + z') 
Jo 



dt 



dz' 



where 



dt\ -1 



!z j = -flo(i + z)yj(i + n mZ )(i + z) 2 - n A (2z + Z 2). (7) 

Hq = lOO/i Km s _1 Mpc -1 is the current Hubble constant. Throughout the paper we adopt a 
flat cosmology with h = 0.65, density parameters fi m = 0.35 and S7a = 0.65, and a Cold Dark 
Matter power spectrum of density fluctuations with power-law index n = 0.96, an amplitude 
as = 0.87 and baryon density parameter Of, = 0.04 (Bahcall et al. 1999). The explosion energy, 
E, is highly uncertain. Based on X-ray afterglow data, Freedman & Waxman (1999) inferred 
explosion energies in the range 10 51 ' 5 — 10 53,5 ergs(A$7/47r), where AO is the solid angle into which 
the energy is channeled. If one assumes that GRB explosions are isotropic, then in some events 
the radiated energy is estimated to be in excess of 10 53 ergs, reaching a value of 3.4 x 10 54 erg for 
GRB 990123 (see, e.g., Kulkarni et al. 1999). It has been suggested that these values are reduced 
by orders of magnitude due to beaming (Kulkarni et al. 1999; Meszaros & Rees 1999), although 
the relatively low efficiency for converting internal shock energy into radiation might imply an 
energy budget as high as 10 54 ergs even after the beaming correction (Kumar 1999). Unless 
otherwise noted, we adopt the value E = 10 52 ergs, which is intermediate for the range inferred by 
Freedman & Waxman (1999). In § 3.3 we will discuss the effect of beaming and different choices 
of the explosion energy on our results. 
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2.2. Non-Relativistic Regime 

As the blast wave decelerates, the fireball eventually makes a transition from relativistic 
to sub-relativistic (sr) expansion (Huang, Dai & Lu 1998; Woods & Loeb 1999; Frail, Waxman 
& Kulkarni 1999; Wei & Lu 1999). In this regime, the evolution of the GRB remnants is well 
described by the Sedov- Taylor solution (Taylor 1950; Sedov 1959). The time when this transition 
takes place in the observer frame, t sr , can be approximately derived by setting the shock velocity 
in the Sedov- Taylor solution to be equal to the speed of light, c. This gives: 



(1 + *) 



E /0.47 



m p n 



1/3 

= 1.8 x 10 2 (1 + z)£^ 3 n7 1/3 days. (8) 



In the sub-relativistic regime, the shock radius and velocity evolve as, 

r(t)=r(t sr )(t/t sr f/\ (9) 

v(t) = c(t/t sr )- 3/5 , (10) 

where all times are in the observer frame and r(t sr ) is the shock radius at time t sr . For a strong 
non-relativistic shock, the post-shock particle density, energy density, and magnetic field are, 

ri = An, (11) 
9 

u = -nm p v , (12) 
8 

B' = (8ire B u') 1/2 . (13) 

The observed flux requires a substitution of the values of u c , v m and F Um in the sub-relativistic 
regime into equations (1) and (2). The synchrotron formulae can be used since the emitting 
electrons are still ultra-relativistic for the relevant remnant ages in this regime. The characteristic 
synchrotron frequency is given by, 

"M^W < 14 > 

Recalling that 7 m oc B' 2 oc u'/n', we find 7 m oc v 2 oc t -6 / 5 , and 7 C oc (B' 2 t)~ l oc i 1 / 5 . Equation (14) 
then provides an expression for v c and v m in the sub-relativistic regime. In particular v c oc £ -1 / 5 
and v m oc t~ 3 . During the relativistic expansion the observed peak flux is constant in time; 
however, in the sub-relativistic regime it varies as F Um oc r 3 B' oc i 3 / 5 . 

In the UV regime one needs to take account of absorption by the intergalactic medium 
(IGM). At redshifts greater than the reionization redshift, z re i on , the neutral IGM is optically 
thick to photon energies above the Lya resonance. Recent models for the reionization of the IGM 
predict z re i on in the range of 7-12 (Haiman & Loeb 1997; Miralda-Escude, Haehnelt & Rees 1998; 
Valageas & Silk 1999; Chiu k Ostriker 1999; Gnedin 1999; Ciardi et al. 2000). Here we adopt 
Zreion = 8; the precise choice of which has only a minor effect on our results since the Lya forest 
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yields a high opacity at the redshifts of interest even if most of the IGM is ionized. For z £ 5, we 
use the updated Lya forest opacity derived by Haardt & Madau (1996) and Madau, Haardt & 
Rees (1999), which is based on the observed absorber distribution in the spectra of high-redshift 
quasars. The continuum depression blue-ward of the Lya resonance is already close to unity at 
z ~ 5 (see Figure 13 in Stern & Spinrad 1999), and so the treatment of the Lya absorption at 
5 ^ z ^ 8 has little effect on our results. 

Next, we derive an expression for the ambient gas density as a function of redshift. For 
simplicity, we assume that the ambient gas has the mean density of the interstellar medium of the 
host galaxy, as inferred for GRB970508 and GRB971214 (Waxman 1997; Wijers & Galama 1999). 

2.3. Evolution of the Gas Density in Galactic Disks 

The phenomenological modeling of many afterglow lightcurves implies values of the ambient 
medium density ~ 1 cm~ 3 , which are of order the mean density of the interstellar medium of disk 
galaxies (Waxman 1997; Wijers & Galama 1999; Meszaros 1999). In popular hierarchical models 
of galaxy formation (Kauffmann, White <fe Guiderdoni 1993; Mo, Mao & White 1998), the mass 
and size of galactic disks evolves with redshift (Barkana & Loeb 2000). Hence, in modeling the 
statistical properties of GRBs within disk galaxies at different redshifts, we need to follow the 
evolution of the average density of the interstellar gas within these galaxies as a function of cosmic 
time. 

We model a disk galaxy as a gaseous disk embedded within a dark matter halo. The disk is 
assumed to be self-gravitating and radially exponential with a scale radius, r d . Its scale height at 
any radius is dictated by the balance between its self-gravity and the gas pressure. For simplicity, 
we assume the disk to be isothermal. Its mass density profile is then (Spitzer 1942), 

p(C,r)=p e- r/rd sech 2 (£/V2), (15) 

where £ and r are the vertical and radial coordinates of the disk, respectively; po =const, 
£ = C(4-7rG / ooe _r ' /rd /Cg) 1//2 , and c s is the effective sound speed (including the possible contribution 
from turbulence). We adopt a sound speed of c s ~ 10 km s _1 , corresponding to a temperature 
of ~ 10 K below which atomic cooling is highly inefficient. At redshifts z ^ 15 the formation of 
galaxies with virial temperatures below 10 4 K is expected to be suppressed due to various feedback 
effects (Barkana & Loeb 1999; Ciardi et al. 2000). 

The scale-radius of the disk is dictated by the angular momentum per unit mass of the gas 
and can be expressed as (Mo, Mao & White 1998), 

r d = -^=^-\r vir . (16) 
V2 m d 

Here j d = J d /J and m d = M d /M, where J and M are the halo angular momentum and 
mass and the subscript d refers to the disk; A and r vir are the halo spin parameter and virial 
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radius respectively. We assume that the disk mass fraction is the cosmic baryonic fraction, 
Md/M = (f^/f2o)> since gas cooling is efficient at the high redshifts of interest. The characteristic 
size distribution of local galactic disks is obtained by adopting jdl m d = 1 (Fall & Efstathiou 
1980). Based on numerical simulations, we parametrize the distribution of the spin parameter A 
by a log-normal form (Mo, Mao & White 1998, and references therein), 



p(\)d\ = — 7n=exp 



In 2 (A/A) 
2*1 



with A = 0.05 and a\ = 0.5. The halo virial radius is given by (Barkana & Loeb 1999), 

M n(z) 200 



r vir — 7.56 

with 5 r = 187T 2 and, 



hr 



1O 8 M n m 5 C 



1/3 

(l + zj^kpc, (18) 



n (z) = ^(i + zf ( i9) 

The total mass of the disk, Md, can be derived by integrating equation (15). This yields the 
relation, 

» " m^ d (20> 

The star formation rate depends on the local surface density of the gaseous disk. Schmidt's law 
implies that the star formation rate is proportional to the local surface density of the disk to the 
power of 1-2 (Kennicutt 1998). Hence we first model the probability for a GRB occurrence within 
a volume element 2nrdrd( as being proportional to the density squared, Pgrb = An 2 ((,r); while 
in §3.3 we consider the sensitivity of our results to other power-laws scalings. The proportionality 
constant A is chosen to normalize the integral of Pgrb to unity. Consequently, we can write the 
average of the ambient density probed by GRBs as 

/•oo /*oo 

(P)= / p((,r)P G RB2Trrdrd( = 0.29p . (21) 
Jo Jo 

The number density of atoms in the medium is thus (n) = (p)/p,m p , where p is the mean molecular 
weight. Similarly, we can derive (^/n) and (1/n) and substitute them into equations (3)-(5). Even 
if the central engine of all GRBs were standard, the dependence of the ambient medium density 
on the properties (mass and size) of the host galaxy would have introduced significant scatter and 
evolution to the flux distribution of GRB afterglows. 



3. Results 
3.1. Flux Evolution 



Figure 1 shows the observed flux for a GRB hosted by the average halo mass as a function 
of redshift. This average is calculated based on the Press & Schechter (1974) formalism, which 
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provides the mass function of dark matter halos as a function of redshift. The two sets of curves 
correspond to two observed wavelengths: in Figure la the thick lines correspond to A Q b s = 2/im 
and the thin lines to 5000A, while in Figure lb the thick lines correspond to A Q b s = 10 cm and 
the thin lines to 1 mm. Within each set, the lines correspond to an observed time of 1 hr (solid), 
1 day (dotted) and 10 days (dashed). While the optical and infrared emission peak at observed 
times shorter than an hour, the radio emission has not reached its maximum even after 50 days. 

At observed frequencies below the Lya resonance frequency, u a (z) = 2.47 x 10 15 /(1 + z) Hz, 
the emitted flux is not absorbed by the intervening intergalactic medium at a redshift z. Figure 1 
shows that at these frequencies the afterglow flux is only weakly dependent on redshift. This can 
be crudely understood from the scaling laws implied by equations (1) and (2), by considering the 
simple Einstein-de Sitter cosmology (£l m = 1), for which oc (1 + z)[l — (1 + z) -1 / 2 ]. Recalling 
that n oc (1 + z) A for a fixed host galaxy mass [see eqs. (18)- (21)] and substituting these scaling 
laws into equations (1) and (2), we find for the case of fast cooling, 

r (i + *) 5 /2[i_(i +z )-i/2]-2 Vc>v , 

F„<xl {\ + z)-*/*[\-(\ + z)-V*\- 2 v m >v>v c , (22) 

[ (l + *)(P-6)/4[l_(l + z )-l/2]-2 v>V m - 

while for slow cooling, 

(l + *)5/6[l_(l + ;8 )-l/2]-2 Vm>Uj 

(l + *)(P-3)/4[l_(l + z )-l/2]-2 Uc>u >V m , (23) 

(1 + z)(P- 6 )/ 4 [l - (1 + z)-V2]-2 V > V C . 

Although these scaling laws ignore the evolution of the characteristic galaxy mass with redshift, 
they demonstrate that the observed flux does not decline rapidly with increasing redshift, but 
might rather increase with redshift at low frequencies (e.g. in the millimeter or radio regimes). 
This behavior results from three causes: (i) for a fixed observed time-lag after the GRB, the 
cosmic time-dilation implies that the higher the source redshift is, the earlier the emission time is, 
and the brighter the intrinsic afterglow luminosity is; (ii) the spectral slope of GRB afterglows 
results in a favorable K-correction for the redshifts and wavelengths of interest; and (iii) the gas 
density in galaxies increases with redshift (due to the corresponding increase in the mean density 
of the universe) . 



3.2. Number Count 

In a snapshot mode of observations, the total number of GRBs with observed flux greater 
than F at an observed frequency u is 

»r, „ x f°° f°° , ,^t F (z,M,v) dn . . dV , % , , 

N(>F,u)= / r GRB (z,M) f \' ' > ( z )—dMdz. 24 

JO JM min (1 + *) dM dz 
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where vgrb denotes the GRB rate in a galaxy halo of mass M at a redshift z, tp is the observed 
time over which a GRB event from this halo yields a flux brighter than F at a frequency v, 
dn/dM is the comoving number density of halos with masses between M and M + dM (based 
on the Press-Schechter formalism), and M m j n is the minimum galaxy mass in which stars form 
at a redshift z. The minimum halo mass inside which star formation occurs, is related to a 
minimum virial temperature of ~ 10 4 K below which atomic cooling of the gas is suppressed 
and fragmentation into stars is inhibited (Haiman, Rees & Loeb 1997; Ciardi et al. 2000). This 
minimum virial temperature implies M m \ n = 4.4 x 10 9 Mq(1 + z) _1 ' 5 /i _1 (Padmanabhan 1993), 
and is also consistent with our choice for the sound speed of the disk (see § 2.3). The comoving 
volume element per unit redshift, dV/dz, is given by, 



dV Aired'- 



dz (1 + z) 



dt 

dz 



(25) 



where dp and dt/dz are given in equations (6) and (7), respectively. Next we consider the value of 
the GRB rate, vgrb- 

If GBRs are related to the final stages in the evolution of massive stars, then one may assume 
that the GRB rate, vgrb, is proportional to the star formation rate (SFR). This assumption is 
justified as long as the time delay between the formation of a massive star and a GRB event is 
short compared to the Hubble time at the redshift of interest. We adopt this as our working 
hypothesis. Wijers et al. (1998) have estimated the constant of proportionality between vgrb and 
the cosmic SFR for an Eistein-de Sitter cosmology, but their results do not change significantly 
for other values of the cosmological parameters (Bagla 1999, private communication). Using the 
observed redshift history of the SFR per comoving volume in the universe (Lilly et al. 1996; 
Madau et al. 1996) and the observed distribution of 7-ray flux of GRBs, Wijers et al. (1998) 
have calibrated the GRB rate per comoving volume. Assuming that GRBs are standard candles, 
they have found a GRB rate at z = of R GRB (0) = 0.14 ± 0.02 Gpc~ 3 yr" 1 . This estimate 
ignores recent corrections to the cosmic SFR. The same calculation, performed with a SFR per 
comoving volume that flattens at high redshift, as suggested by more recent corrections for dust 
extinction, still leads to similar results (Bagla 1999, private communication). Assuming that the 
proportionality constant does not vary with redshift, we write, 

r GRB (z,M) = RGR f { ^ M^z,M). (26) 

Here p*(0) is the star formation rates per comoving volume in units of M yr^ 1 Mpc~ 3 , and 
M+(z,M) = a(z)r2{,M is the SFR within a particular halo in units of M© yr" 1 . We calibrate the 
star formation efficiency as a function of redshift, a(z), so as to match the observed cosmic SFR, 
p*, at z 4, 

P*(z) = / M*(z) — (z)dM. (27) 
JM min dM 

An analytical fit for the observed p*(z) (assuming a Salpeter IMF and an extinction correction of 
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Ai500 = 1-2 mag) was derived by Madau & Pozzetti (1999), 

2Se 3 4z 

^) = e 3. 8z + 44 , 7 Moyr^Mpc- 3 . (28) 

At z > 4 the SFR within a particular halo is calculated using the common prescription 
a(z) oc tjy n (z), where td yn (z) is the mean dynamical time of the halo. The constant of 
proportionality is chosen so as to match the observed cosmic SFR at z = 4 (see Ciardi et al. 2000, 
for more details). 

Figure 2 presents the total number counts based on equation (24). The set of curves 
corresponds to different observed wavelengths. From the right to the left, A j, s is equal to 10 cm, 
1 mm, 2 fj,m and 5000 A. At low fluxes the number counts approach an asymptotic value. Out 
of the entire population of faint GRB afterglows, we predict that there are ~ 15 GRBs from 
redshifts z £ 5 across the sky which are brighter than ~ 100 nJy at an observed wavelength of 
~ 2(um. Figure 3 illustrates the contribution to the total number counts from different redshift 
bins, centered at z=2 [Fig. (3a)], 6 [Fig. (3b)] and 10 [Fig. (3c)]. Infrared or radio afterglows 
are detectable out to redshifts as high as 10, while optical afterglows are strongly absorbed by 
the intervening IGM at z ft 5. Figure 4 shows the differential number count distribution per 
logarithmic flux interval and redshift interval of all sources at a given redshift. At high frequencies, 
for which the observed flux decreases with increasing redshift, low-redshift events dominate the 
counts at most fluxes. However, at low frequencies (i.e. millimetric wavelength) the high-redshift 
events dominate the counts at sufficiently low fluxes. 



3.3. Dependence of Results on Model Assumptions 

Next we consider the dependence of our results on various model parameters. Figure 5a 
shows the observed flux at A {, s = 2^m and t = 1 day for different choices of the mean ambient 
density. The standard case with Pgrb oc n 2 (solid line) is compared to two other cases; one 
where Pgrb oc n (dot-dashed line) and the other where the ambient density is kept constant 
and equal to 1 cm~ 3 (dotted line). The figure indicates that for observations in the infrared, all 
three cases provide identical results at z £ 5, when one enters a regime where the observed flux 
becomes independent of the ambient density [either the fast cooling regime with u > u m or the 
slow cooling regime with v > v c , given by Eqs. (1) and (2)]. Outside this regime the observed 
flux is roughly proportional to n 5 / 6 and yjn in the fast and slow cooling cases, respectively. In 
the slow cooling regime, for example, the condition v > v c leads to a minimum redshift above 
which the flux is independent of ambient density (1 + z) > 0.32i?^ 1 nj" 2 t^ 1 (A b s /2/um) 2 . Hence, at 
long wavelengths or very low densities, the density independent regime is not reached. Figure 5a 
shows that the difference between the Pgrb oc n 2 and Pgrb oc n cases, is small since they follow 
a similar evolution of the gas density with redshift. Figure 5b shows the number count of GRBs 
per logarithmic redshift interval at z = 6, both for the standard case with Pgrb oc n 2 (solid line) 
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and the case of a constant ambient density equal to 1 cm 3 (dotted line). The difference between 
these cases is small. 

Figure 6a shows the dependence on the energy output of the afterglow number count from a 
redshift bin centered at z = 6 for an observed wavelength of 2/um. We consider two cases that 
bracket the standard case (E = 10 52 ergs); namely E = 10 53 erg (dotted) and E = 5 x 10 50 erg 
(dashed). The increase in the observed flux as the energy output increases [see eqs. (l)-(5)], is 
reflected in the afterglow number count. 

So far we have assumed that the GRB explosions are spherically-symmetric. If, however, the 
energy release is beamed, then the afterglow lightcurve is expected to evolve differently than we 
assumed. At early times, as long as the expansion Lorentz factor, 7, is still larger than the inverse 
of the jet opening angle, 6* _1 , the expansion behaves as if it were spherically symmetric with 
the same energy output per solid angle. Once the jet decelerates to a Lorentz factor, 7 £ 1 , 
the lightcurve declines more rapidly, since the jet starts to expand sideways and to reduce the 
mean energy output per solid angle (Rhoads 1999a,b; Panaitescu & Meszaros 1999). Finally, the 
isotropization of the energy ends when the expansion becomes sub-relativistic, at which time the 
lightcurve recovers the spherically-symmetric Taylor-Sedov evolution for the actual total energy 
output. (Synchrotron radiation losses are usually negligible in the late afterglow phase.) For an 
observing point which is aligned with the GRB jet axis, the afterglow lightcurve may start with 
an isotropic-equivalent energy of E and end in the non-relativistic regime with the actual energy 
output of E sr = 7]E where 77 = (Afi/47r) = 2(ir9 2 /An) < 1 is the fraction of sky around the GRB 
source which is illuminated by two opposing jets of angular radius, 9. Hence, the afterglow makes 
a transition between the lightcurve corresponding to a highly-relativistic explosion with a total 
energy E to the lightcurve of a non-relativistic explosion with an energy r\E. To examine the 
effect of beaming on our results, we model this transition using a linear interpolation between 
the two lightcurves. We start with an isotropic equivalent of E = 10 53 ergs, and adopt 6 = 0.1 
so that the actual energy output is E sr = 5 x 10 50 erg (since 7 ^ 10 2 S> initially). As 
expected, the afterglow fluxes and source number counts make a transition between their values 
for the two boundary energies. Figures 6a-c show the transition between the two regimes with a 
dot-dashed line. We find that, for A Q bs = 2/xm, the number counts are dominated by sources in 
the relativistic phase. However, at longer wavelengths, X b s = 1 mm or 10 cm [Fig. 6b or Fig. 6c, 
respectively], the bulk of the emission takes place at late times when the expansion is in its mildly 
or sub-relativistic phase. Obviously, the effect of beaming can be recovered from the shape of 
the afterglow lightcurve in individual GRB events. Figure 6 shows that the existence of beaming 
might also be inferred from the number count statistics. 

4. Discussion 

We have calculated the expected fluxes and number counts of high-redshift GRB afterglows, 
assuming that their rate of occurrence is proportional to the star formation rate. We have 
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computed the observed properties of GRB lightcurves at different wavelengths ranging from the 
optical to the radio, treating both the relativistic and sub-relativistic phases of the expansion of 
their remnants. 

Our main result is that at a fixed observed time after the GRB event, the characteristic 
afterglow flux is not decreasing rapidly with increasing redshift (Figs, la, lb), in contrast with 
other high-redshift sources, such as galaxies or quasars. Hence, the broad-band spectrum of GRB 
afterglows is ideally suited for probing the ionization state and metal content of the intergalactic 
medium at high redshifts, particularly during the epoch of reionization (Loeb 1999; Lamb & 
Reichart 1999). 

The only difficulty in using GRBs as probes of the high-redshift universe is that they are rare, 
and hence their detection requires surveys which cover a wide area of the sky (see the vertical axis 
in Figures 2-4). The simplest strategy for identifying high-redshift afterglows is through all-sky 
surveys in the 7-ray or X-ray regimes. In particular, detection of high redshift sources will become 
feasible with the high trigger rate provided by the forthcoming Swift satellite, to be launched 
in 2003. Swift is expected to observe and localize ~300 GRBs per year, and to repoint within 
20-70 seconds its on-board x-ray and UV-optical instrumentation for continued afterglow studies. 
The high-resolution GRB coordinates obtained by Swift, will be transmitted to Earth within 
~50 seconds. Deep follow-up observations will then become feasible from the ground or using 
the highly-sensitive infrared instruments on board the Next Generation Space Telescope (NGST), 
scheduled for launch in 2008. Swift is sufficiently sensitive to trigger on the 7-ray emission from 
GRBs at redshifts z £ 10 (Lamb & Reichart 1999). We note that even if some GRBs occur outside 
of galaxies, the intergalactic medium would be sufficiently dense to produce an afterglow at these 
high-redshifts since its mean density is £ 10~ 4 [(1 + z)/10] 3 cm -3 . 

The effect of beaming is expected to be minor for sources with millimeter and radio fluxes 
^ 10~ 5 Jy, since the number counts at these fluxes are dominated by sources which are in the 
sub-relativistic phase of their expansion, and for which the energy output in the explosion is 
already isotropized (see Fig. 6b, 6c). 

For a characteristic energy output of GRBs of ~ 10 52 ergs, our model implies that at any time 
there should be ~ 15 GRBs from redshifts z ^ 5 across the sky, which are brighter than 100 nJy 
at an observed wavelength of 2 fim. The Next Generation Space Telescope (NGST) will be able to 
measure the spectra of these sources. Prior to reionization, the spectrum of GRB afterglows might 
reveal the long sought-after Gunn-Peterson trough (Gunn & Peterson 1965), due to absorption by 
the neutral intergalactic medium. 

We thank F. Haardt for providing us the IGM opacity curves, and J. Bagla and R. Croft for 
useful discussions. BC thanks the CfA pre-doctoral fellowship program for support during the 
course of this work. This work was also supported in part by NASA grants NAG 5-7039 and NAG 
5-7768, NSF grant 9900077, and by a grant from the Israel-US BSF. 
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Fig. 1. — (a) Observed flux for a GRB hosted by a "typical" halo (having the mean mass) as 
a function of redshift (see text for details). The curves are for a frequency v = 6 x 10 14 Hz 
(^obs = 5000 A, thin lines) and v = 1.5 x 10 14 Hz (A & s = 2^m, thick lines). From the top to the 
bottom different observed times after the GRB are shown: one hour (solid line), 1 day (dotted) 
and 10 days (dashed). 
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Fig. 1. — (b) As in Fig. (la), but for v = 3 x 10 11 Hz (X bs = 1 mm, thin lines) and v = 3 x 10 9 
Hz (X bs = 10 cm, thick lines). The dot-dashed line corresponds to an observed time of 50 days. 
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Fig. 2. — Total number of GRBs with observed flux greater than a limiting flux, F, at different 
observed wavelength \ bs- From the right to the left X Q b s is equal to 10 cm, 1 mm, 2 fim and 5000 

A. 
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Fig. 3. — (a) Contribution to the total number count (Fig. 2) of GRBs from a redshift bin centered 
at z=2. The observed wavelengths are, from the right to the left: X b s = 10 cm, 1 mm, 2 /im and 
5000 A. 
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Fig. 3. — (b) As in Fig. (3a) with a redshift bin centered at z=6. 
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Fig. 3. — (c) As in Fig. (3a) with a redshift bin centered at z=10. 
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Fig. 4. — (a) Number count distribution per logarithmic flux and redshift interval. The observed 
wavelengths are, from the right to the left: \ Q b s = 10 cm, 1 mm, 2 /im and 5000 A. The redshift 
bin is centered at z=2. 
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Fig. 4. — (b) As in Fig. (4a) with a redshift bin centered at z=6. 




Fig. 4. — (c) As in Fig. (4a) with a redshift bin centered at z=10. 
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Fig. 5. — (a) Observed flux for a GRB hosted by a "typical" halo (having the mean mass) as a 
function of redshift (see text for details). The flux is shown for the reference observed wavelength 
2 //m and an observation time of 1 day after the GRB. The different curves correspond to different 
ambient densities, namely: a constant density of 1 cm -3 (dotted line), and density profiles of galactic 
disks with Pgrb °c n 2 (solid line; standard reference case) and with Pgrb oc n (dot-dashed line). 
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Fig. 5. — (b) Contribution to the total number count of GRBs from a redshift bin centered at z = 6 
for an observed wavelength of 2 fim. Notation is the same as in Fig. (5a). 



-27- 



10° E 



10 



-1 _ 



CV2 
I 

-a lo- 2 



10" 





N 




+ 


A 










o 




9 



10 



-4 _ 



10" 5 



i i i nnij — i i i inii| — i i i niiij — i i i iiiiij — i i i iiiiij — i i i iinij — i rrmg 

(a) Khs= 2 M n : 

z = 6 



E 



53 



E 53 & = 0.1 



J 50.7 \ 



\ 



I I I I I I i i 1 1 



10" 8 10" 7 10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 
Limiting Flux [Jy] 



Fig. 6. — (a) Contribution to the total number count of GRBs from a redshift bin centered at z = 6 
for an observed wavelength of 2 /im. The different curves refer to different energy outputs: E = 10 53 
erg (labeled -E53; dotted line), and 5 x 10 50 erg = 10 50 ' 7 erg (labeled £50.7; dashed line). Two of the 
curves refer to a beamed explosion with an equivalent isotropic energy output of 10 53 ergs and a 
beaming angular radius of = 0.1 (dot-dashed line; see text for details). The actual energy output 
in this case is 5 x 10 50 erg. 
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Fig. 6. — (b) As in Fig. (6a), but for the reference observed wavelength of 1 mm. 
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Fig. 6. — (c) As in Fig. (6a), but for the reference observed wavelength of 10 cm. 



